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Some Nonnegative Trigonometric Polynomials Connected 

With a Problem in Probability 

Eugene Lukacs^ and Otto Szdsz^ 

Let 0<^6i<C62<C . . . <C^n be integers and let g{6) be the Vandermonde determinant 
formed from 6i, 6i, . . . , 6^ with the first row replaced by sin^ {bid/2) (i=l, . . . , n). The 
function g(0) is then a cosine polynomial. In connection with a problem in probability, the 
question arose as to when g(6) is a nonnegative trigonometric polynomial. This question is 
answered for six classes of such trigonometric polynomials. 



In a previous paper ^ we gave a necessary condi- 
tion that a polynomial without multiple roots must 
satisfy in order that its reciprocal be the Fourier 
transform of a distribution function. Imposing a 
further restriction on the polynomials it is easy to 
derive the following condition: 

The reciprocal of a polynomial whose roots are 
all single and have the same imaginary part is the 
Fourier transform of a distribution function if and 
only if 

(1) The polynomial has one purely hnaginary root 
ai (a 7^0) and n pairs of com.plex roots ± b^+ai 
(0<b,<b2< . . . <bn, k = \, 2, . . ., n). 

(2) The determinant 
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>0 for all 19. 



This condition follows easily from fonnula (4.4) of 
the reference cited in footnote 3. It is therefore 
of some interest to study tliis determinant and to 
investigate for what values of bi, 62, • . ., &« it is 
a nonnegative function of 6. 

In tliis paper we consider this determinant only 
for integer values of the bt and show that it repre- 
sents for certain configurations of the bi a nonnegative 
trigonometric polynomial. Certain relations for 
generalized Vandermonde determinants of odd in- 
tegers are also obtained. 

We introduce first some notations. Let 
0<6i<62<C . . . <^bn be n integers and 
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the generaUzed Vandermonde determinant formed 
from the numbers bl, b^ . . ., bl \\\i\\ the exponents 
k, 1, 2, . . ., {n-l). Clearly ¥^ = ¥2= . . . = 
Vn-\ = 0, while Fo is the ordinary Vandermonde 
determinant of the numbers b\, . . ., bl. Denote 
further by A,„ the minor of the element in the first 

n 

row and m-th cohunn of T so that ^0= S (~ 1 Y~^^m • 
In this paper we study the function 
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For the discussion of g(d) we need the following 
lemmas: 

Lemma 1, If 61, 69, . . ., 6n are integers, the 
determinant (2) can be factored so that 



/7((9) = (1— cos ey yl(cos <9), 



(3) 



where Ai^x) is a polynomial in x of degree bn—n. 

Lemma 2. If 61, 60, . . ., ^n are odd integers we 
have 

.^'(i9) = sin(9(l-cos2^)^-i7i(cos^) = (sin(9)2"-7?(cos(9), 

(4) 

where B{x) is a polynomial in x of degree 6,^—2^^ + 1. 

To prove lemma 1 we differentiate (2) with respect 
to d and then set ^=0. This shows that 






for j=l, 2, ... ad inf 
fori=0, 1, 2, . . . (71-1) W5) 
fort=7?, (71+1), • . . ad inf. ^ 
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From (5) we obtain the expansion of g{d) into a series 

(6) 



9{e)=j:i-if-' '" 



k=n 






On the other hand, (2) indicates that g{d) is a cosine 
polynomial of degree b„ which vanishes for 6=0. 
Therefore it has the form 



^(0)=^i(l-cos e)+^2(l-cos dY+ 

+ ;8.(1-COS0)*» 



(7) 



Comparing (6) and (7) it is seen that pi= ^2= . . . = 
i8„_i=0 so that (1 — cos 6)^ is a factor of g{d); this 
establishes lemma 1. 

To prove lemma 2 we assume that all the bi(i=ly 
2j , . . J n) are odd numbers. It follows then from 
(2) that 

g(T-e)=2Vo-g{d), (8) 

therefore 

g'{T-d)=g'{d), (9) 

From (3) we have 

g'{d)={l — cos ey-'^ sin e[nA(cos d) — 
(1 — cos (9) A' (cos d)]. 

Substituting this into (9) we obtain 

(1 — cos ey-'^ sin ^[71^ (cos ^) — (1 — cos ^)^'(cos 6)] 

= (l + cos ey-^ sin einA(—cos 6) — 
(l + cosl9M'(-cos^)]. 

This shows that 7i^(cos ^) — (1 — cos e)A^(cos 6) and 
therefore also g^{e) has the factor (1 + cos eY~^. This 
completes the proof of lemma 2. 
If we introduce as a new variable 



and write 

P(x)=g{eiTG cos x) 
we have from lemma 1 

P{x) = {l-xyA(x), 
and similarly from lemma 2 

P'(x) = -il-x'y-'B(x), 



(10) 



(11) 



(12) 
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We substitute next the expansion (6) into (8) and 
see that 



^"•-S'-^-OT! 



6»2*= 
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(14) 



Differentiating (14) and setting ^=0 we obtain from 
this equation a number of relations for generalized 
Vandermonde determinants formed of odd integers 



2Fo: 



m T/ 

^„^ '> (2k)\ 
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for m=l, 2, 



(2n-l) 



fc=7i+m 



S (- 
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(2A:-27i-2m)! 

for m=:0, 1, 2, . . . ad inf. 
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for m==0, 1, 2, . 



= 



ad inf. ^ 



In the following we discuss several configurations 
of the integers 61, 62, . • ., ^n^ which lead to non- 
negative trigonometric polynomials g(Q), the results 
are given in statements, labelled (A) (B) . . . (F). 

(A) If the hi are the first n consecutive integers, that 
is \ihi^=i for i = l, 2, . . ., ti, then the trigonometric 
polynomial g{Q) is nonnegative. 

Prooj: From lemma 1 we see g{d)=A{l —cos dY, 
where ^ is a constant. Hence ^=2~^^(7r); from (2) 

n 

it is seen that ^(x) = Fo+X)^m. Therefore A>0 
and consequently g{d) nonnegative for all values of 6. 

(B) If the hi are the first n consecutive odd integers, 
that is, if hi=2i—l for i = l, 2, . . . , ?^ then the 
trigonometric polynomial g(d) is nonnegative. 

Prooj: Since 6^=27i--l the polynomial B{x) of 
lemma 2 reduces to a constant B. From (13) we have 
P\x)=-B{l-x^Y-'' or, in view of P(l)-0, 



P{x)=B 



Cii-ty-' 



-'dt. 



Formulae (11) and (2) show then that P(0)=^(7r/2) = 
Fo^O, hence 



B= 



P{0) 



j\i-tr 



->o. 



'dt 



Therefore P(x)>0 for |x|<l. This completes the 
proof of statement (B) since the inequalities 
P{x) >0 for \x\ < 1 and gld) >0 for all d are equivalent. 

(C) If hn=n-\-l, that is if the numbers 61, . . .,6^ 
are obtained from the first consecutive (n+ 1) integers 
by omitting the integer k (l<k<n) then the trigo- 
nometric polynomial g(d) is nonnegative if and only 
if 2P>(7i+l).. 

Proof: In this case P(x) has degree (71+ 1) so that 
according to lemma 1 A(x) is a linear function. 



A(x)=a-{-hx. 



(16) 
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The substitution (10) ti'ansforms cos kd into Ttc{x) 
where Tu{x) is the k-ih Tchebichcff polynomial of 
the first kind. Formulae (2), (3) and (16) show that 

P{x) = {l-xY{a-^hx) = V,+ ±,{-ir^mT,^{x). (17) 

m = l 

The coefficients a and h of A{x) can be determined 
from (17), after some elementary computations it is 
seen that 

^W = 2''A„[.+^±l)±^]. (18) 

In the interval \x\<^l the functions A{x) and P(x) = 
(1 —xy A{x) have the same sign so that it is sufficient 
to determine when A(x) is nonnegative in |x|<Cl. 
From (18) it is seen that A(x) >0 in fx|<l if and only 
if n(n+l)+k^>{n+iy -k^, that is if 2F>n+l. 

(D) If the bi (i=l, 2, . . ., n) are all odd numbers 
and if bn=2n-\-l so that 6i, 62, . • -, bn are obtained 
from the ffi-st (n+1) odd integers by omitting the 
k-ih. odd integer 2k — I (l<k<n), than the trigono- 
metric polynomial g{d) is nonnegative if either one 
or the other of the following conditions is satisfied 



but 



k>^(l + ^/T+2n) 



t<^(i+vr+2^), 



(i) 
(ii) 



n — 2^(^—1) 



2[n(n+l)-k{k-l)] 
where Zn is the root of the equation 



<2^ 



/-' 



{l-tY-\f-z')dt=0 



(19) 



which falls into the interval (0, 1). 



If neither (i) nor (ii) is satisfied, then the function 
g(d) assumes also negative values. Moreover it is 
possible to sim.plify (ii) for large n by proving the 
Corollary to statement (D). 

If for large n and k<C.[l + -\ll+2n]/2 also 



V' 



(n+i)[n-2^(^-l)] 



<P = 



=0.5939157. 



2[nin+l)-k{k-l)\ 
then g(6) is nonnegative; if on the other hand 



(n+\yn-2kik-l)] 
2[n{n+l)-k(k-l)] 



>P' 



then there are integers n and k for which g{d) assumes 
positive and negative values. 



To prove statement (D) we need the following 
lemma 

Lemma 3. Let 



and 



Ps{x)= r(l-tV-H''dt for 5 = 0, 1 (20) 



Q{x,n=P,(x)-t'Po(x). (21) 



The polynomial Q(x, f) is nonnegative in the interval 
— 1 <x< 1 if and only if 

f<^n. (22) 

Here Zn is the root of the equation 
R(2)^Q{-z,z)^ r (l-ty-\t'-z')dt=0, (23) 

wliich falls into the interval (0, 1). 

Moreover let Wn^-^n^—Zn- 
Then lim Wn = p exists and is the root of the equation 

n—*oo 

F(p)--=V^(^l-p'^-2p'f\l-t')e-f>'t^dt=0, (24) 

which is located in the interval f Oj-W^ )• 

Some values of Zn as well as p were computed. 
The proof of this lemma is rather lengthy. In 
order to avoid interrupting the discussion of the 
various nonnegative trigonometric polynomials, this 
proof will be given in the last section of this paper. 

We proceed now to the proof of statement (D). 

According to lemma 2 B{x) is a polynomial of the 
second degree. Since the bi are all odd numbers 
P^{x) and B(x) are even functions of x so that one 
obtains from (13) and (2) 

P\x) = -{l-xy-\a + cx') = ^{-l)'^AmT:jx), 



This relation permits us to determine a and c, and it 
is seen that 



P'{x) = 
-2'''{2n+l)An{l-xy-' 



n-2k{k-l) 



2[n(n+l)-k(k-l)] 
If we use for brevity the notation 

y=2'''(2n+l)An 



and 



_ n-2k{k-l) 
^"•' 2[n{n+l)-k{k-l)]' 



(25) 
(26) 
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we have 



P{x)^y ^\l-t'Y-\t'-tU)dt\ 



(27) 



If (i) is satisfied we see that f,^,A;<0, so that from (27) 
P'{x)<{) for |a:|<l. In this case P{x) is non- 
increasing in the interval — 1 <x<-\-\, since P(1) = 
this means P(x)>0 for |:r|<l so that also ^(^)>0. 

We assume next that (ii) is satisfied, then ff,,yt>0. 
Since 7>0 and P(x)=yQ{x, ^n,k), we see that state- 
ment (D) follows immediateh^ from lemma 3. 

(E) If the numbers bi, . . . , bn are obtained from 
the first consecutive (n-j-2) integers by omitting the 
two integers k and p where l<k<p—l<n, then 
g{d) is nonnegative if and only if a certain polynomial 
Q(x) of the second degree is nonnegative for |x|<l. 
Here 

Q{x) = 4:x'+4:[7i-L{n, k, p)]x+2n''+{n-2)- 

4.nL(n, k, p)+N{n, k, p), (28) 



where 



L{n,k,p)-- 



_{n + 2)[(:n + iy-k'][{n + iy-p']') 



[{n+2y-k'][{n + 2y-p'[ 



in + 2){2n+3){n^-k^){n'- p') 
^^Vh'<^,P) [(n + 2y-k'][{n+2y-p'] 



y (29) 



Proof: In this case bn—n-\-2, so that according to 
lemma 1 A{x) is a polynomial of the second degree, 
A(x)=ax^+bx+c. From (12), (2) and (3) we have 



P{x) = (l — xY{ax'+bx + c) 

n 

= Vo+j:(-l)'"AmT,Jx). 

m = l "" 



(30) 



From this relation we can determine the coefficients 
a, b, and c. A somewhat tedious elementary 
computation shows that 

A(x) = (ax'+bx+c) = 2''-''AnQ(x), 



This completes the proof of statement (E) which 
follows immediately from (30). 

Certain particular cases can be discussed easily 
{£i)lik=p—l=n then N(n, n,n+l) = L(n, n,n-{-l) 
= so that Q(x)=4:x'+4:nx+2n^+n-2. The dis- 
criminant of Q{x) is then —16(n — l){n-{-2) and is 
negative for n.>l. Therefore, Q(x)>Q for all x so 
that g(d) must be nonnegative. 

(b) li k<Cn while p = (n+l) then L(n, k,n + l) = 
and 

Q{x) and its discriminant can be computed easily 
and it can be shown that g{d) is nonnegative if and 
only if 

j^,^ {n-l)in+2) 

(c) li k=n — l and p=n then A^(/?,7/ — /,/? ) = and 

_nin + 2) 



L{n, n — 1, n)^ 



3(n+l) 



The discriminant of Q(x) is then 



16 



9(n+l)= 



(-2n^-57iH7/i2 + 277i+18), 



which is negative if n>3 so that in this case ^(^) >0. 
If 71 = 2, Q{x) has two roots in the interval 
— I<x<+1 so that g{d) assumes also negative 
values 

We proceed to discuss a more complicated case 
by assuming that the 6i, b2, . . . , b^ are obtained 
from the first (n + 2) odd integers by omitting two 
odd integers {2k —1) and (2p —I) where 



l<k<p-l<n 



(31) 



so that 6„=2«, + 3. 

We see then from lemma 2 that B{x) is a poly- 
nomial of degree four. By the procedure employed 
in case (C) it is possible to determine the coeffi- 
cients of this polynomial. 



We obtain 



Here 



P'(x)=-y{l-x-)''-'B{x)=-y{l-xy-\4x' + ax--b) 



(32) 



7 = 2^"(27i + 3)A„>0, 

_ 2 {-2in+iyi3n+5) + {n+l){2n+5)[kik~l)+p{p-l}]~4kik-l)pip-l) 
" [{n + 2)in + l)-k{k-imn + 2){n+l)-pip-l)] 

{~(n+l)(3n + 4) + 2{n+l)[kik-l)+p{p-l)]-4k{k-l)p{p-l)} 
[{n + 2){n+l)-k{k-l)][in + 2){n+l)-p{p-l)] 



(33) 
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The trigonometric polynomial g{6) is nonnegative if 
and only if 

P{x)=^y r il-tT-\4t'+at'-b)dt>0 for !il\\x\<l. 

(34) 

The nonnegativity of the polynomial P{x) depends 
on the nature of the roots of the polynomial 
BixJ—Ax'^-i-ax^—b. This introduces the need to 
distinguish a number of cases. For this purpose 
we write 

A(r) = 4T^+ar — 6 = 4(r — Ti)(r — 72) where \ 



166 f 



-a—^a'+lQb _-a+^Ja'+16b ?' (^^) 
^'~ 8 ^'- 8 ) 

We are now in a position to formulate statement (F) . 
(F) If the numbers 61, 62, • • •, ^n are obtained 
from the first consecutive (n+2) odd integers by 
omitting the two integers 2fc— 1 and 2p—l (where 
l<k<p — l<n) then the trigonometric polyno- 
mial g{6) is noimegative if and only if one of the 
following four, mutually exclusive, conditions is 
satisfied 

(a) a>0 andb<0, 

(0) a<0 anda2+166<0, 

(7) «^^0 and 0<J<4+a and also 

i?2(V^)-TA(V^)>0. 

(d) 0>a> — 4 and — --<6<a+4 and also 
Id 

R2(^^r'2)^rJil(^2)>0, and if ti>0 

/?2(-A/rO-T27?i(-V^)>0. 

Here 

R,{b)= r {l-ty-H'^'-'\f-z')dtioTi=l,2, (36) 

so that Ri{z)=R{2) as given by (23). 

Proof; If A(r)=4T2+aT-6>0 for 0<r<l, then 
P\x)^ -y{l-x'Y-^h{x')<^ for -l<x<l. Since 
P(1) = we see that P{x)>{) for -l<:r<+l. 

Clearly /i(t)=( 2r+T) "~T^""^j therefore 

hir) > in < r < 1 if and only if (37) 

(^2r+fy>^^±AMinO<r<l. (38) 

If a>0, (38) holds if and only if 6<0, but this is 
exactly condition (a) . 



We next derive the condition for the validity of 
(37) when a<0. Clearly A'(r)-=8T+a, A''(r)=8, so 
that h{T) has a minimum at T=—a/S. 

We consider first the case |a|/8= —a/S<^l, we 
have then 

mm h{T) = h {^—J=- 



0<r<l 



16 



Relation (37) is therefore satisfied if —S<ia<iO and 
a'+lQb<0. If |a|/8>lweseethatA'(r)<8(r-l)<0 
in 0<r<l, so that min A(T)=A(l)=4+a— 6. 

0<r<l 

The function /i(r) is nonnegative if a< — 8 and 
4: + a-b>0. Buta2+166=(a+8)2-16(4+a-6)for 
any a, hence if a^+ 166 < Owe always have 4 4-a—6>0. 
The alternatives discussed for case a<^0 can there- 
fore be joined so that condition (/3) is established. 
To complete the discussion we must consider the 
remaining cases, namely 



and 

(50 



a>0 and 6>0, 

a<0 anda2+166>0. 

In case (7') we see from (35) that ti<0 while t2>0. 
Since P (x) =47 r (1 -tY~' if - n) {f - T2) dt, it fol- 
lows that P(x) <0 if T2> 1, that is if ^la'-\-l^b>a+S, 
but this occurs if and only if 6>4+a. 

Thus P{x) <0 for -l<x<l ifa>0 and 6>4+a, 
therefore (7') reduces to the case 



(7") 



a>0 and 0<6<4+a. 



In this case both roots ri and T2 are real and a 
simple computation shows that min P{x)=P( — V'^2), 

hence P(x)>0 in |x|<l if and only if 
That is if 



f- 






TiC ^{l-tT-'{t'-r2)dt>0 

J--\T2 

or, written in a more concise notation 

SO that (7) is established. 

We proceed with the discussion of the last case 
{d^)a<CO and a^+166>0 and show first that always 
a+4>0. From (33) we obtain easily 



(a + 4) ^2{ n+iy(n'+n-l) + {n-\-l)[k(k-l)+p(p-l)]-2k{k-l)p(ip-l) 
2 [^n + 2){n+l)-k{k-l)][{n+2){n+l)-p{p-l)]^ 
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On account of (31) this fraction has always a pos- 
itive denominator so that its sign is determined by 
the sign of the numerator which we denote by 
M{n,k,p) . We have then 

M{n,k,p)=p{p-l)[(n+l)-2k(k-l)] + 

{n+l)k{k-l) + 2{n+iy(n'+n-l) 

= k{k-l)[{n+l)-2p{p-l)] + 

{n+ l)p(p)- 1) + 2(71+ 1)2(^2+^- 1). 

If (n+l)-2k(k-l)>0 then M(n,k,p)yO. If on 
the other hand (n+1) —2k{k —1)<C0 then also 
(n+l) —2p{p—l)<i0. Thereiore M(n,k ,p) decreases 
for fixed k and increasing p as well as for fixed p and 
increasing k. Hence M{n,k,p) attains its minimum 
value if k and p are as large as possible, that is 

M{n,k,p) >M(n,n,n+l) 

= 2{n+iy{n^+n-l) -2n\n+l){n -2) 

>2n\n+iy -2n\n+l)(n -2) ii n>l. 

Therefore M{n,k,p)>6n\n+l)yo if n>l so that 
also a+4>0. It is therefore no restriction if we 
write the case (5') in the form 

(8') 0>a>-4anda2+166>0. 

We must first consider the case 

(di) 0>a>— 4 and b >a+4. 

We have then a^+166>0 so that h(r) has two real 
roots; from (35) we see that ti<^0 while t2>1. 
Therefore 

P{x) = ^y r {l-ty-\t'-Ti){t'-T^dt<0 

for — 1 <a::<l so that (8i) does not yield nonnegative 
trigonometric polynomials. We finally have to in- 
vestigate the possibility that 



(« 



0>a>-4 and-Y6<&<4 + a 



holds. Again a^+166>0 and A(t) has two real 
roots Ti and t2. It follows then from (35) that 
l>T2>ri>-l. 

Assume first that ri>0. 

By a simple computation it is seen that P{x) has 
the two minima+VrT and— Vrl inside the interval 
(— 1,+ 1). P(x) is therefore nonnegative for 
-1 <a;< 1 if and only if P{^f7^>0 and P(-V^) >0. 
According to (34) and (35) this means that 



J- 



-V^ 



{l-tT~'t\t'-T2)dt' 






and 



J -v^ 



ri)dt- 



T2 r {l-tT~\t'-Ti)dt>0 



or, using the notation of (36) 

i?2(V^-T,J?i(V^>0 

and simultaneously R2{.—^fr^— r2Ri{—^fu) ^0. 

If ri<0 then there is only one minimum for P{x)j 
and therefore only the first condition remains. 
This establishes condition (5) and completes the 
proof of statement (F). 

In this section we give proof of lemma 3, which 
was stated on page 141. We write 



P,{x)^ Cil-tV-H'^dt {s = 0, 1) 



(^9) 



and want to determine the conditions which f has 
to satisfy in order that the polynomial 



Q{x,^)=P,(x)-t'Po(x) 



(40) 



should be nonnegative for — l<a:;< + l. 

Clearly 0<P^{x)<Po(x) for -l<a:< + l. There- 
fore Q{x,^)<CQ for |a;|<l if f >1. If for some value 
fo such that |fo|<Cl the function Q(Xj^q)>0 for 
|x|<l then also Q{x,t)>0 for |a;|<l and |f!<|fo|. 

To obtain a criterion for the nonnegativity of 
Q(Xy^) in \x\<l we have to determine the greatest 
possible f such that |f|<l and Q{x,^)>0 for \x\<l. 

An elementary computation shows that the func- 
tion Q(x,^) has exactly one minimum in the interval 
— I<a:;< + 1 which is located at x=— f. 

We consider next the equation R{z) = Q(—ZjZ)=0 
or, written in greater detail, 

r {l-tT-\t'-z'^dt==Pi(-z)-z^Po{-z) = 0. (41) 

Clearly i?(0)>0 while i?(l)=Pi(-l)-Po(-l)<0 so 
that the equation (41) has at least one root in (0,1). 
Moreover 



'^=-2^j\i-tr-^dt 



so that dR/dz vanishes in the interval 0<2<1 only 
at the point 0=0, therefore the equation (41) has 
exactly one root Zn in the interval (0,1). 

Since min Q(XjZn)=R{Zn) = we see that Q{x,Zr) 

-l<a:<l 

>0 in |a:|<l and therefore Q(a?, f)>0 in la;|<l 
if f <^n. 

If however f >^w then 

e(-0,, ^)=P,{-Zn)-i'P,{-Zn) 

= -{^^-zl)Po{-Zn)<0 
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so that Q{Xj f) assumes also negative values. This 
establishes the first part of lemma 3. We still have 
to discuss the asymptotic behaviour of the solution 
of equation (41). 

We derive first a useful inequality. By definition 
Zn is a root of equation (41), therefore 



If we set 

then we see that 
2 



Wn = -Jn + - Zn 



(42) 






di>{). 



or 



We see therefore that 

rwr(l) ^. r(.)r(l) 
r(.+|) n+^r(n+i) 



>o 



so that 



y.<^orzl<^- 



(43) 



In order to derive an asymptotic estimate for Wn 
we rewrite (41) by dividing the interval of integra- 
tion into two parts (—2,0) and (0,1). The integrals 
over the second interval may be easily expressed in 
terms of gamma functions. A simple computation 
yields : 



If we write 



2(n+|)r(.+i) 



i(2) = 



Biz) 



and 



Tin) 
h{z)=ril-2Hr-\l-t')dt, 



(44) 



we have 



Liz)=^^^[l-i2n+l)z']- 



2r 



^(n+'-)zU 
Tin) y'+2)^^'^- 



From the asymptotic formula for log r(x) it is seen 
that 

{n+-yT{n) 



hence 

L{z)=^ll-i2n+l)2^]-2eo( 



(45) 



We next transform the integral /«. Let 0<x<l 
then (l-x)6Xl and (l-x)^-^<e-^'^-^)^ therefore 

0<l-[{l-x)eT~'<{n-l)[l-{l-x)e']. 

Since e''>l + x we have also (1— a:)e^>l — x^ and 
l-(l-a:)6^<x2 so that finally 0<l-(l-x)^-^e^^-^>* 
<i{n-l)x^ or 0<e-^''-'^^-{l-xY-'<(n-l)x\ We 
replace here x by zH^ and see that 0<e-(«-i)2^^^— 

(l-32^2)n-l_g_(^_l),2^2_^^_-^)^^^4^4f()j. ^2^2^1 

with 0<6>X1. 

We substitute this into the expression (44) for In 

and obtain 

I^=\ e-(n-i)z^t^(l-t')dt-^{n-l)z' (46) 
Jo ^^ 

for 22<i with 0<5X1. 
We have 

-^22,2 3 , / — 5-22,2 \ 

Since 0<l-6 ^^ <|^2^2 ^e see that (l-^ ^ J 
= rjzH^ with 0<77<3/2 and hence 

therefore 

C\-(n-i)z^t^(l -t')dt= r e~^'''^^^'''\l -V)dt 
Jo Jo 

+ ^ zh-in-Dz^. 
15 



with 0<o'<l and 0<X<-- 
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Substituting this into (46) we have 



2X^%_,,_,,,._25.(n-l)3^ 



15 



with 0<8n<l, 0<o-<l, 0<X< 



35 



(47) 



3 



We introduce next a new variable Wn 



--■^n + ^z, 



using (45) and (47) we see that Wn satisfies the 
equation 



■VTT 



-2wl 6^(^-') 



2d^wt (n-l) 



2x wi ^ n+y 

' — • e 2 

15 , 1 



35 



= 



(48) 



For each value of n this equation has a root Wnj by 
(43) we see that the set of these roots is bounded. 
We consider any accumulation point p of this set 
and a subsequence [Wn^ of the sequence {Wn} such 
that lim Wn^^p- From (48) it is seen that p satisfies 

the equation 

i^(p)=V^(^-p^)-2p3|^-p^^^ {\-mt=0 (49) 

and by transforming the integral we obtain finally 
for p the equation 

F{p)=^(~-p'^-^\p'-x^e-^>dx = Q (50) 

To complete the proof of lemma 3 we have only to 
show that the sequence { Wn ] converges. We dem- 



onstrate this by showing that [Wn] has only one 
accumulation point, and this follows if we show that 

F(p)=0 has only one root 



in(o,^i)_ 



We see easily that J^(0)>0 and ^(a/^)<0 so 

that there is at least one root in the interval. 
From (50) we see that 

F\p) = -2p[^JT+2 re-''dx'] 

so that F\p)<0 for p>0. This shows that F{p) 
has exactly one root. 

The equation (41) has been evaluated for certain 
values of n, similarly the asymptotic solution p has 
been found from (50). This was done in part at the 
Computation Laboratory of the National Bureau of 
Standards, Washington, D. C, and in part at the 
computation department of the Institute for Numeri- 
cal Analysis in Los Angeles, Calif. The results are 
given in the following table, which shows also the 
values of Wn. 



n 


Zn 


Wn 


2 


. 38197 


. 60394 


3 


. 32114 


. 60079 


4 


. 28244 


. 59915 


5 


. 25505 


. 59814 


6 


. 23434 


. 59745 


7 


. 21798 


. 59696 


8 


. 20463 


. 59659 


9 


. 19346 


. 59630 


10 


. 18395 


. 59606 


15 


. 15122 


. 59536 


20 


. 13141 


. 59500 


p = 0.5939157 



Washington, May 25, 1951 
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